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ABSTRACT 

We investigate tlie effect of a global change in the ionizing continuum level on the behavior 
of the strong optical broad emission lines seen in spectra of the nuclear emission-line regions of 
active galactic nuclei (AGN), including the Balmer lines, He I A5876, and He II A4686. Unlike 
most of the prominent heavy element lines found in the UV, the optical hydrogen and helium 
recombination lines' emissivities are strongly dependent on the incident continuum flux, since 
these lines arise out of excited states whose optical depths depend on the incident flux of pho- 
tons. Using photoionization calculations we determine the luminosity-dependent responsivities, 
r]{r, Lit)) = A log L;,;„e /A log Lcont, of these lines for a general model of the broad emission line 
region (BLR), with the purpose of establishing them as important probes of the physical con- 
ditions within the BLR of AGNs. The dependence of these lines' emissivities on the incident 
photon flux invokes a dependence in their responsivities on distance from the central continuum 
source. In particular, the responsivities of these lines are generally anticorrelated with the in- 
cident photon flux. Thus, their responsivities vary with distance within the BLR for a fixed 
continuum luminosity and change with time as the continuum source varies. Consequently, after 
correcting for light-travel-time effects the response of the Balmer and optical helium lines should 
generally be strongest during low continuum luminosity states. Responsivity that depends on 
photon flux and continuum state may explain a number of outstanding problems currently under 
investigation in broad-line variability studies of these and other emission lines. These include 
the origin of the intrinsic Baldwin effect, measurements of luminosity-dependent lags (a "breath- 
ing" BLR) and luminosity-dependent variations in: integrated broad emission-line flux ratios 
(including He II A4686/H/3), broad line profile shapes, and radial velocity-dependent intensity 
ratios. The broad Ha/H/? and H/3/He I flux ratios and the Balmer emission-line responsivity are 
observed to decrease from the line center to the line wings. These, along with our findings, lead 
to the conclusion that the BLR velocity field diminishes with increasing distance from the central 
continuum source. This is consistent with recent reverberation studies that find a relationship 
between the emission-line lag and rms profile width for multiple lines in individual AGN, which 
implies that the velocity field is dominated by a central massive object. Finally, the responsivity 
of ionization-bounded clouds can account for much of the observed behavior of the optical re- 
combination lines (e.g., the weak response of the Balmer line wings) previously attributed to a 
substantial contribution from matter-bounded clouds at small BLR radii. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert — (galaxies:) quasars: emission 
lines — line: formation — line: profiles 



1. INTRODUCTION 

Historically, attempts at interpreting the ob- 
served variability behavior of the broad emis- 
sion lines in active galactic nuclei (AGN) have 
tended to regard the broad emission line region 
(BLR) as a stationary entity (non-evolving). By 
this we mean that any variability in the observed 
emission-line intensities has been attributed solely 
to reverberation (light-travel-time) effects within 
a spatially extended BLR. Thus, variations in the 
measured BLR sizes and the amplitude of the 
line response have been assumed to arise primar- 
ily through differences in the temporal variability 
of the ionizing continuum, which manifests as a 
difference in shape and FWHM of the continuum 
autocorrelation function (ACF) from one event to 
the next (e.g., Perez, Robinson, & de la Fuenta 
1992). In essence, higher frequency continuum 
variations are better able to probe smaller BLR 
sizes, while larger BLR sizes are better probed by 
lower frequency continuum variations. 

While a difference in the continuum variability 
time scale (e.g., the width of the continuum ACF) 
explains in part some of the observed variability 
seen in the broad emission lines in AGNs, it is by 
no means the complete picture. Indeed, more so- 
phisticated photoionization calculations (e.g., see 
Krolik et al. 1991; O'Brien et al. 1995; Bottorff 
et al. 1997; Kaspi & Netzer 1999; Horne, Korista, 
& Goad 2003) have modeled the line variability 
by including the local emission-line response to 
changes in the overall luminosity and/or shape of 
the ionizing continuum. The models presented by 
Carroll (1985) and Netzer (1991) anticipated the 
non-linear response of some of the emission lines. 

In this paper we illustrate, using a general de- 
scription of the BLR, how large changes in the 
mean continuum level may explain much of the 
observed variability behavior of the optical re- 
combination lines. In § 2 we describe the grid 
of photoionization models used in this study. In 
§ 3 we describe the physical origin of line rcspon- 
sivity and its impact on the following issues re- 
lated to BLR variability: the origin and slope of 
the intrinsic Baldwin effect (§ 3.2; e.g., Kinney, 
Rivolo, & Koratkar 1990); variations in the mea- 
sured continuum-emission-line time-delays (§ 3.3; 
e.g., Peterson et al. 2002), variations in the emis- 
sion line flux ratios (§ 3.4; e.g., Balmer decrement: 



Tran, Ostorbrock, & Martcl 1992; H/3/He II 4686: 
Peterson & Ferland 1986); changes in proflle shape 
(§ 3.5; e.g.. Wanders & Peterson 1996); and 
velocity-dependent line intensity ratios (§ 3.6; e.g., 
Stirpe, de Bruyn, & van Groningen 1988). We 
emphasize that the variations in these quantities 
discussed here are not a consequence of reverber- 
ation effects within a finite-sized BLR, although 
they are intimately related. Rather, they are due 
to global changes in the mean ionization state of 
the BLR gas. In § 4 we discuss the impact of these 
findings on the purported contribution of a very 
broad optically thin component to the emission- 
line response, and other related issues. It is our 
hope that many of these observed emission-line 
variations and trends reported in the literature 
may be examined in a new light. We summarize 
our main findings in § 5. 

2. PHOTOIONIZATION MODELS 
2.1. The Grid of Photoionization Models 

For the purposes of demonstration, we adopted 
the grid of photoionization computations pre- 
sented in Korista & Goad (2000; hereafter KGOO) 
that they used to model the strong UV broad 
emission lines of NGC 5548. The reader is di- 
rected there for details; here we mention some 
of the salient features of their grid and adopted 
model. 

Using Ferland 's spectral synthesis code, Cloudy 
(v90.04; Ferland 1997; Ferland et al. 1998), KGOO 
generated a grid of 3249 photoionization mod- 
els of broad line-emitting entities, here assumed 
to be simple, single total hydrogen column den- 
sity slabs {Nh = 10^^ cm~^), each of which 
has constant gas density (for simplicity) and a 
clear view to the source of ionizing photons. 
The grid dimensions spanned 7 orders of mag- 
nitude in total hydrogen gas number density, 
7 < logn// (cm~"^) < 14, and hydrogen-ionizing 
photon flux, 17 < log$i/(cm~^ s~^) < 24, and 
stepped in 0.125 decade intervals in each dimen- 
sion. We refer to the plane defined by these two 
parameters as the density-fiux plane. The choice 
of spectral energy distribution (SED) is described 
in KGOO (see also Figure Al in Horne, Korista, 
& Goad 2003). In the Appendix we discuss the 
sensitivity of the results presented here to these 
and other assumptions. 
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In the top four panels of Figure 1 we show the 
equivalent width (EW) distributions of the promi- 
nent optical broad emission lines of Ha A6563, H/? 
A4861, He I A5876, and He II A4686, namely, those 
of primary interest in the present paper. These 
EWs are all referenced to the same incident con- 
timmm flux at I2I5 A and are a measure of the ef- 
ficiency by which ionizing continuum photons are 
converted into line photons. The peak in the EW 
distribution for each line is denoted by a triangle. 
Note that the EWs of the Balmer and He I lines 
increase with increasing density, especially at the 
lower values of the photon flux, illustrating the ef- 
fects of coUisional excitation from excited states 
(Ferland & Netzer 1979). In every case, the line 
EWs decline in the direction of increasing ionizing 
photon flux for values greater than those marking 
the peak EWs: log^jj = 17.125, 18.250, 19.125, 
and 20.750 for Ha, H/9, He I, and He II, respec- 
tively. The rapid decline in the line EWs of the 
two Balmer lines for log^jj » lognH-|-10.25 is due 
to the hydrogen in the fixed-column density slabs 
becoming fully ionized, a subject we discuss later. 
The two helium lines do likewise at slightly larger 
values of log(UHc) = log$^f — lognn- The upper 
left corners of these diagrams represent gas near 
the Compton temperature. Detailed discussions of 
these and other EW contour plots in the density- 
flux plane can be found in Korista et al. (1997). 
In order to aid in comparison and to serve as a 
fixed point of reference, the star symbol in each 
of the panels shows the location of the standard 
BLR parameters from Davidson & Netzer (1979). 

The EW distributions of two other hydrogen 
recombination lines, H7 A4340 and Lya A1216^, 
are plotted in the bottom panels of Figure 1. We 
refer to them from time to time for purposes of 
comparison. 

2.2. A Simple Model for the BLR of 

NGC 5548: A Test Case 

To model the mean UV spectrum of NGC 5548 
from the 1993 FST campaign (Korista et al. 1995; 
hereafter K95) KGOO took the "locally optimally- 
emitting clouds" (LOG) approach of Baldwin et al. 
(1995) by summing the emission from a weighted 
distribution of "clouds" along the two dimensions 

This is a sum of Lyo and minor contributions from O V] 
A1218 and He II A1216. 



of gas density and radius; KGOO assumed con- 
stant column density clouds of 10^'^ cm~^. The 
cloud distribution functions in gas density and ra- 
dius are described in KGOO. In brief, the model- 
integrated spectra included clouds spherically dis- 
tributed around the continuum source whose gas 
densities^ and distances from the central source of 
ionizing photons spanned 8 < lognn (cm~^) < 
12 and 1-140 It-days; very highly ionized clouds 
{logUHC > 11.25) were excluded from the integra- 
tion (see KGOO for details). Asymmetric geometric 
effects, such as emission-line anisotropy (Ferland 
et al. 1992; O'Brien et al. 1994), were not consid- 
ered. KGOO showed that this simple model with 
widely distributed cloud properties was also gen- 
erally successful in reproducing the observed UV 
broad emission-line variability. 

With the choice of cosmological parameters 
and adopted Galactic reddening along the line 
of sight to NGC 5548 from KGOO, the mean value 
of AFa(I350) continuum flux from the 1993 HST 
campaign (K95) corresponds to log ALaisso (ergs s~-^) 
« 43.54. Using the central continuum SED 
adopted by KGOO, the corresponding mean hydrogen- 
ionizing luminosity from the 1993 HST campaign 
is log Lion (ergs s^^) « 44.26. At this luminosity a 
hydrogen-ionizing photon flux log <f>//(cm^^ s^^) = 
20 corresponds to a distance from the continuum 
source of i? ~ 12.6(75/i?o) It-days, and the model 
BLR spans log$ff(cm-2 s'^) w 17.9 to « 22.2 in 
this mean continuum state. In this paper we con- 
sider two historical extrema in the AFx(1350) con- 
tinuum flux of NGC 5548 about the above mean 
level. In 1992, the continuum went into a very 
low state, log ALai350 (ergs s~^) ~ 42.90, while a 
high state of log ALai350 (ergs s^^) « 43.81 was 
recorded during the 1989 campaign. This is a fac- 
tor of about 8.2 change in the continuum flux at 
1350 A, and we assumed that the SED did not 
change shape. 

While we do present results based on a simple 
LOG model of the BLR in NGC 5548, they should 
be representative of most general classes of models 
describing the broad emission lines. The major 
trends should be largely model- independent, while 
the details should be diagnostic of the physical 

^ Given the choice of a single cloud column density and 
the assumption of constant hydrogen gas densities within 
the cloud, the cloud radial thicknesses range over 10^^ — 
10^5 cm. 
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conditions within individual BLRs. 

3. EMISSION-LINE RESPONSIVITY 

Summing the hnc emission along the cloud gas 
density distribution at each radius, we can derive 
an effective cloud surface flux of the emission lines 
at that radius (see KGOO; Baldwin et al. 1995). 
Figure 2 shows the surface flux of H/3 as a func- 
tion of radius from the central continuum source 
for each of the two continuum states. The slope 
of cloud emission-line smface flux plotted against 
the ionizing flux (oc ijo„/r^) is a measure of the 
effective radial responsivity of the emitting gas for 
small amplitude (or equivalently, short time scale) 
changes in the ionizing continuum (see Goad et al. 
1993). Similarly, the difference in the two surface 
flux curves plotted in Figure 2 (in logj^g space) rel- 
ative to the logarithm of the factor of 8.2 in the 
incident continuum flux measures a "global" re- 
sponse to long-term, large amplitude, changes in 
the continuum state. We define this responsivity 
parameter as 



V{r) 



Alog F{r) 



line 



Alog$H 



(1) 



This definition of line responsivity converges to the 
definition given in Goad et al. (1993) for very small 
variations in the continuum flux^. Note that the 
upper curve in Figure 2 is essentially the lower 
curve shifted along the x-axis by 0.5 x log(8.2) 
decades. A logarithmic slope of —2 in Figure 2 
corresponds to ri{r) = 1, whereas flatter slopes 
correspond to ri{r) < 1. When integrated over the 
full line-emitting region, what we call the effective 
responsivity, r/e//, is the logarithmic slope in the 
relation Lii„e oc L^^^t ■ 

The solid curves in Figure 3 illustrate the global 
responsivities as functions of radius of H/3 (bottom 
left), Ha {top left), He I A5876 {top right), and 
He II A4686 (bottom right) for this large change 
in continuum state. We emphasize that the major 
changes in responsivity with radius are not a result 
of our choice of fixed-column density clouds; this 
is discussed below. The dashed and dotted curves 
represent the responsivities for small changes in 



*In an alternative approach to analyzing the emission-line 
response to continuum variations, Sparke (1993) compared 
the width of the emission-line ACF with that of the con- 
tinuum ACF. 



the continuum flux about the high and low contin- 
uum states, respectively, and these are discussed 
in § 3.1 and § 3.3. It is clear from the global re- 
sponsivities that larger relative variations of these 
emission lines' surface fluxes due to changes in 
the continuum state can be expected to occur at 
larger radii where the incident continuum fluxes 
are smaller. 

We can understand the general trends in line re- 
sponsivity, with a minimum number of model de- 
pendent assumptions, through a study of Figure 1. 
We can relate the responsivity to the gradient in 
the EW contours in Figure 1 as 



dlogEW 
dlog^H 



1. 



(2) 



The responsivity is proportional to the efficiency 
in converting a change in the number of ionizing 
continuum photons into escaping line photons, and 
so is in proportion to the change in EW contours 
shown in Figure 1 relative to the incident photon 
flux. For a fixed scale in cloud distance from the 
central ionizing source, the EW contours in Fig- 
ure 1 shift vertically upward with a decrease in the 
continuum luminosity, and shift vertically down- 
ward with an increase in the continuum luminos- 
ity. Clouds in the density-flux plane for which the 
line EW distribution is either flat or changes little 
with the incident photon flux have responsivities 
that lie near 1 (i.e., roughly 1:1 linexontinuum 
variations) for that line (see, especially. He II 
A4686). Those clouds whose line EWs increase 
with increasing incident photon flux have respon- 
sivities exceeding 1 for that line (bottoms of the 
panels for each of these lines except Lya). Those 
clouds whose line EWs decrease with increasing 
incident photon flux have responsivities less than 
1 for that line: the steeper the gradient, the lower 
the responsivity. Finally, where the line EWs de- 
cline very rapidly with incident photon flux in the 
density-flux plane (near the central diagonal of the 
figures), the emission line responsivities become 
negative. 

To illustrate these effects, in Figure 4 we super- 
pose a responsivity grey-scale on the EW contours 
for each of the lines in Figure 1. The grey-scale 
representations are: 77 < (white), < < 0.5 
(light gray), 0.5 < r] < 1.0 (medium gray), 1.0 < 
ri < 1.5 (medium dark gray), 1.5 < r/ < 2 (dark 
gray), and r/ > 2 (black, appears only in He II). 
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Given that we generally measure < rief / < 1 for 
the broad emission lines, this diagram shows how 
the emissivity and responsivity together constrain 
the cloud properties, especially their distribution 
in radius from the continuum source. 

We now briefly describe how the responsivities 
of the optical recombination lines arise. At suf- 
ficiently high gas densities and ionizing photon 
fluxes, each of these four emission lines' EWs di- 
minish with increasing value of the ionizing photon 
flux (Fig. 1). This indicates a lessening efficiency 
of the conversion of ionizing continuum photons 
to escaping line photons for those clouds lying 
nearer to the continuum source. This is primarily 
due to the increasing populations and thus optical 
depths in the excited states of hydrogen and he- 
lium with increasing continuum flux (Netzer 1975; 
Netzer 1978; Ferland & Netzer 1979; Ferland, Net- 
zer, & Shields 1979; Kwan 1984; Rees, Netzer, 
& Ferland 1989; Ferland et al. 1992; Shields & 
Ferland 1993), i.e., thermalization at large opti- 
cal depth. At a given continuum flux, the opti- 
cal depth of Ha is largest, followed by H/3, He I 
A5876, and He II A4686 (Rees et al. 1989). Various 
background opacities also play a role (Shields & 
Ferland 1993), and photoionization out of excited 
states of hydrogen becomes effective for densities 
exceeding ~ 10^" cm""^. In fact, those high den- 
sity clouds lying above roug hly the EW(Ha) = 5 A 
contour (0.7 dex) in Figure 1 lack a hydrogen ion- 
ization front for this reason. As a consequence 
of these optical depth effects, the more substan- 
tial relative flux variations in these four emission 
lines tend to be weighted by clouds lying at larger 
radii. This latter point should be true to some 
extent for most emission lines, since they too suf- 
fer from diminished efficiency in the conversion 
of ionizing continuum photons into escaping line 
photons at smaller radii where the gas densities of 
clouds that emit the lines must be higher. Refer 
to Figure 2 of KGOO for the EW contours in the 
density-flux plane for the stronger UV emission 
lines. The Mg II A2800 emission line is most like 
the Balmer lines in this respect, and we would ex- 
pect that line to have a similarly low responsivity 
(see Figure 2b in KGOO). This is indeed generally 
observed (e.g., see Clavel et al. 1991). Clouds fully 
ionized in hydrogen (such as those lying near or 
within the rapidly declining EW contours running 
diagonally in the density-flux planes in Figure 1), 



with their very low (^ 0) or even negative respon- 
sivities, might in principle be arranged in such a 
way as to substantially reduce the line responsivity 
at the larger radii. However, these low efficiency 
clouds will not contribute much light and there- 
fore will add little to the integrated line response 
unless their covering fractions are very large. 

For time-scales of the order of the light-travel 
time across the BLR, the power spectrum of the 
continuum fluctuations coupled with light-travel 
time effects across a finite geometry will also fac- 
tor into the variations of an emission line's flux 
and perhaps profile;. However, it is the respon- 
sivity that establishes the outer envelope of these 
variations. Light travel time effects can only act to 
reduce or dilute the observed line response. We re- 
iterate that the variations in the line intensity ra- 
tios and line profiles presented here are not due to 
light travel time effects (e.g., see Robinson, Perez, 
& Binette 1990; Perez, Robinson, & de la Fuente 
1992; Pogge & Peterson 1992) or other geomet- 
rical effects (e.g., line anisotropy, O'Brien, Goad, 
& Gondhalekar 1994), but are rather due to dif- 
ferences in the time-averaged photon flux inci- 
dent on a spatially extended BLR and their ef- 
fect on the line emissivity. Krolik et al. (1991), 
Pogge & Peterson (1992), and Gilbert fc Peter- 
son (2003) have found that reverberation effects 
on the Liine oc L^^nt relation can be largely re- 
moved by correcting the emission-line light-curves 
for their mean time-delays. If this is indeed the 
case, then any residual line - continuum flux cor- 
relations must be primarily due to the responsiv- 
ity of the gas in the BLR, as we have deflned here. 
However, if a line's transfer function has a long 
tail in time delay (in comparison to the important 
continuum variability time scales), simply remov- 
ing the lag may not sufficiently correct for this geo- 
metric dilution of the line response. The observed 
lag-corrected emission-line response will then be a 
function of the line's responsivity, as defined here, 
and the prior continuum history. More advanced 
methods, such as that introduced by Hornc et al. 
(2003), may then be needed to disentangle the ef- 
fects. 

Finally, it is apparent from Figs. 1,3, and 4 that 
He II A4686 should be the most responsive emis- 
sion line of the four optical emission lines discussed 
here, followed by He I, H/3, and lastly Ha. This 
is observed (e.g., see Wamsteker et al. 1990; Diet- 
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rich ct al. 1993; KoUatschny 2003). The model el- 
fective global responsivities of these four emission 
lines are listed in column (2) of Table 1. To illus- 
trate the trends within the Balmer line series, we 
also plot the EW distribution of H7 A4340 within 
the density- flux plane in Figure 1, list its effective 
responsivity in Table 1, and illustrate its respon- 
sivity in the density- flux plane in Figure 4. Lya 
A1216 is also a recombination line, but with dif- 
ferent physics, and so we provide the same infor- 
mation for it. 

3.1. The Emission-Line Response Versus 
the Local Ionizing Photon Flux and 
Continuum Luminosity State 

Figure 3 shows that the responsivity of H/3 is 
expected to be a strong function of radius, varying 
from ~ 0.2 in the inner BLR to ~ 1 at the outer 
radius. More generally, this is indicated by the 
gradient in the EW contours with respect to the 
ionizing photon flux in Figure 1 (see also Eq. 2), 
and coarsely illustrated in Figure 4. Another in- 
teresting point raised by Figure 3 is that the H/3 
line response should be significantly stronger dur- 
ing low continuum flux states than during high 
ones for similar continuum fluctuation amplitudes. 
This is illustrated in Figure 3 by the dotted and 
dashed curves representing the low and high con- 
tinuum state H/? line responsivity, respectively, 
computed for small (0.1 dex) variations about the 
high and low states. The effective responsivities 
in the high and low states are given in columns 
(3) and (4) of Table 1. K. Home (2002, private 
communication) and Cackett & Horne (2005) re- 
covered the luminosity-dependent delay map for 
the 13 yr monitoring campaign of the H/3 broad 
emission line of NGC 5548, and found a similar re- 
lation between the amplitude of the line response 
with continuum state. 

Figure 3 and columns (3) and (4) of Table 1 
shows that all of these lines display generally larger 
effective responsivities at larger distances (lower 
incident photon fluxes) and during lower contin- 
uum states. Additionally, the amplitude of the 
change in effective responsivity with continuum 
state diflFers substantially between the Balmer and 
helium lines, and this should be diagnostic of the 
physical conditions within the BLR. By contrast, 
Lya shows very little variation in responsivity 
with cloud parameters (Fig. 4), and as a conse- 



quence shows a very modest change in rjeff with 
continuum state (Table 1). 

We elucidate the change in the radial respon- 
sivity with continuum state shown in Figure 3, 
using He II A4686 by way of example. In both 
continuum states, the responsivity at large radii is 
nearly independent of radius and continuum state. 
This is simply a reflection of the general lack of 
dependence of the line EW with incident contin- 
uum flux for ^ 20.5 (see Figure 4). This is 
in contrast to the case of the Balmer lines, espe- 
cially Ha. Note also in Figure 4 that most of the 
variations of 77(He II) with incident photon flux 
occur for ^ 20.5; here in particular r] declines 
with increasing ^h- This dependence of the re- 
sponsivity on incident continuum flux makes this 
line's radial responsivity at small distances quite 
dependent on the continuum state as well. These 
effects are also illustrated in Figure 3. The "knee" 
in the radial responsivity of He II occurs at the 
radius corresponding to ~ 20.5 for each con- 
tinuum state, or logr w 16.25 (7 It-days) in the 
mean continuum state. Simply put, strong anti- 
correlations between the line EW and the incident 
photon flux will correspond to anti-corrolations 
between r] and the incident photon flux, which 
in turn will result in similar anti-correlations be- 
tween 77 and the continuum state. Likewise, rjeff 
will anti-correlate with the continuum state, as 
weighted by the line's radial luminosity function. 
Based on single power-law cloud distribution mod- 
els, O'Brien, Goad, & Gondhalekar (1995) also 
found "a decrease in the mean response when the 
continuum flux increases." 

By definition, a time-varying responsivity 
r](r, Lcontit)) implies that the line transfer func- 
tion, ^{t,v), is time-dependent (i.e., non-stationary). 
In this case the time dependence is due to changes 
in the time-averaged photon flux incident on a 
spatially extended BLR. We therefore expect the 
measured lags of the emission lines and their pro- 
file shapes to depend on the luminosity state of 
the continuum (§ 3.3 & § 3.4). We note that 
unless the velocity field is self-similar in radius 
(Goad et al. 1999), a profile shape that depends 
on the continuum luminosity state fails a hystere- 
sis test (see Wanders 1994; Perry, van Groningen, 
& Wanders 1994). This is true without invoking 
a redistribution of material within the BLR or 
other geometrical effects as suggested by Wanders 
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& Peterson (1996). 

3.2. The Intrinsic Baldwin Effect 

Detailed analyses of the AGN broad emission- 
line response to continuum variations in spcc- 
trophotometric monitoring campaigns show that 
in all cases the responsivity in most measured 
lines is < 1, even after accounting for light travel- 
time effects (e.g., Krolik et al. 1991). Referring to 
Eq. 2, we see that this requires the line equiv- 
alent width to diminish with increasing contin- 
uum flux and vice versa: the intrinsic Baldwin 
effect {EWune ^conv Kinney, Rivolo, & 
Koratkar 1990: KroHk et al. 1991; Poggc & Pe- 
terson 1992; Peterson 1997; Gilbert & Peterson 
2003). In concurrence with Gilbert & Peterson 
the present work suggests that the origin of the 
intrinsic Baldwin effect is due either in whole or 
in large part to the line's responsivity being less 
than 1 throughout the majority of a distribution 
of largely ionization-bounded clouds (as shown in 
Figure 1), i.e., without the necessity of a large pop- 
ulation of matter-bounded clouds (Shields, Fcr- 
land, & Peterson 1995), although we do not rule 
these out. As discussed above we expect r]eff < 1 
for most lines, and thus most lines should demon- 
strate an intrinsic Baldwin Effect, since the intrin- 
sic Baldwin effect index /3 is related to the respon- 
sivity as /? = rjeff — 1. Based on single power 
law cloud distribution models, O'Brien, Goad, & 
Gondhalekar (1995) showed that an intrinsic Bald- 
win effect can be induced in any line by changing 
its responsivity with continuum level. While we 
concur with these findings, we find a more gen- 
eral and significant origin for this effect, namely 
rjeff < 1 as defined and described in § 3. 

Our simple model for the BLR of NGC 5548 
predicts an rjefj for H/3 of 0.54 in the 1989 high 
state and 0.77 in the 1992 low state (columns (2) 
- (4), Table 1). That is, we expect rjeff and so 
the slope in the intrinsic Baldwin effect to be a 
function of continuum state. Compare the dotted 
and dashed lines in Figure 3. This predicted range 
in rjeff approximately brackets that measured in 
the same object by Gilbert & Peterson (2003) of 
0.53 — 0.65.^ A line's equivalent width is usually 

^Thc logarithmic slope 0.65 was obtained when the value of 
the background galactic continuum at 5100 A was fixed to 
the observed value (Romanishin et al. 1995), while 0.53 was 



measured with respect to the local continuum, and 
Gilbert & Peterson (2003) found that the underly- 
ing continuum becomes bluer when brighter in this 
object, such that ^5100 oc Lf-^^Q. They then sug- 
gested that since the A1350 continuum lies closer 
to the ionizing continuum, it might be better to 
reference the H/3 flux to it, and found the relation 
L{Hj3) oc -£'1350"°'''^^. Our simulations assumed 
that the SED remains constant as its amplitude 
changes. Thus, while it is difficult to make de- 
tailed comparisons, we expect that SED changes 
are a secondary effect on the line responsivity; we 
discuss the impact of a variable SED on line re- 
sponsivity in § 4.2. The point remains that line 
responsivity, as defined here, is likely to be respon- 
sible for much of the observed intrinsic Baldwin 
Effect. 

What can the intrinsic Baldwin effect tell us 
about the BLR? Becatise the EW contours of the 
Balmer lines in particular in Figure^ 1 are so nearly 
proportional to the incident continuum flux, we 
suggest that measurements of the responsivities 
of these emission lines will point to the distribu- 
tion of the optical recombination line-emitting en- 
tities in the density-flux plane. Low density gas 
{uh < 10"'^'^ cm~^) illuminated by low values of 
the ionizing photon flux ($// < 10^® cm^^s^^), 
and gas at any density illuminated by even lower 
values of (< 10^^'^) will emit these optical 
recombination lines, but they will not respond ap- 
propriately to variations in the ionizing continuum. 
The responsivities of these clouds in these lines are 
~ 1 and > 1, respectively. Such models would fall 
far short of reproducing the observed variations in 
H/3, independent of any considerations of the line 
luminosity or lag. Next, the observed mean re- 
sponsivity of the optical He II A4686 fine, rj ^ 0.9 
(see, for example, Dietrich et al. 1993; also Kro- 
lik et al. 1991), necessarily implies that a signifi- 
cant fraction of this line forms in gas exposed to 
high incident ionizing continuum fluxes on average 
(logoff > 20.5 cm~^ s~\ or i? < 7 fight days in 
NGC 5548), as well as relatively high densities^. 
This is just as its measured lag indicates (Clavel 
et al. 1991; K95). Thus, Figures 3 and 4 predict 

obtained when this background light was allowed to vary 
in the fit. 

Higher column density clouds would allow for somewhat 
lower gas densities under this condition, and vice-versa 
(e.g. , see Korista et al. 1997) . 
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a positive correlation between the lag of each of 
these lines and their responsivity, 7?e//, for a par- 
ticular continuum state (see also Carroll 1985). So 
together the luminosity, lag, and responsivity (in- 
trinsic Baldwin Effect) of each of the broad emis- 
sion lines for a range of continuum states should 
constrain the distribution of BLR cloud proper- 
ties in gas density, distance from the ionizing con- 
tinuum source, and column density. The quasar 
tomography method of Horne et al. (2003) is par- 
ticularly well suited to this task, and we discuss 
this further in § 4.3. 

3.3. The "Breathing" Broad Line Region 

Netzer & Maoz (1990) reported that differ- 
ent lags for the broad emission lines of Lya and 
C IV were found for each of the three contin- 
uum variability events during the 1989 campaign 
of NGC 5548 (Clavel et al. 1991). Based on data 
from an earlier optical campaign (Netzer et al. 
1990), Netzer & Maoz also reported a significantly 
different lag for H/3 than did Peterson et al. (1991). 
While acknowledging that a variable continuum 
ACF will change the lag, Netzer & Maoz (1990) 
suggested that at least part of the origin of these 
time-dependent lags lay in a non-linear emission- 
line response. Over the full 13 yr monitoring cam- 
paign of NGC 5548 Peterson et al. (2002; see also 
Peterson et al. 1999) found that the H/3 emission- 
line lag scales with the mean continuum level, such 
that r oc F^{qq, albeit with much scatter. This 
positively correlated behavior of the lag with con- 
tinuum state has been called "breathing." Af- 
ter positing that the UV continuum is more in- 
timately connected to the driving ionizing con- 
tinuum and accounting for the relation that they 
found between the A5100 and A1350 continuum 
bands (i^5ioo « -^1350)1 the relat ion T oc Fjj'y^ re- 
sulted. Peterson et al. (2002) then proposed that 
T oc -Fcont raight be a natural consequence of a 
simple ionization parameter-type argument that 
r oc Cackett & Horne (2005), too, have 

found evidence of breathing in their transfer func- 
tion analysis of the 13 yr spectroscopic dataset of 
Peterson et al. (2002). A breathing BLR has also 
been proposed by Goad et al. (1999) to explain 
line profile variation differences between C IV and 
Mg II in NGC 3516. In this subsection, we address 
a likely origin of breathing BLRs and comment on 
the conclusions drawn by Peterson et al. (2002) in 



regards to this issue. 

The single power law cloud distribution models 
of O'Brien, Goad, & Gondhalekar (1995) predicted 
a relation between line lag and incident contin- 
uum level (for a fixed SED) due to a non-linear 
emission- line response. The results presented here 
for the optical recombination lines based on more 
general models arrive at the same conclusion. The 
shift in the lag as a function of the continuum 
state is a reflection of an outward shift in the 
line's luminosity-weighted radius to gas at larger 
radii that is better able to reprocess the increased 
continuum luminosity into line emission, and vice 
versa in low continuum states. The EW contours 
in Figs. 1 and 4 shift vertically as the source lu- 
minosity changes, causing corresponding changes 
in the radial responsivity and hence responsivity- 
weighted radii. We therefore expect that elevated 
continuum states will be accompanied by an in- 
crease in the line lag and a decrease in the line 
responsivity. 

Breathing will not occur for lines for which 
the responsivity as a function of radius and ion- 
izing continuum level is constant. A constant 
ri{r, L{t)) implies that the emissivity-weighted and 
responsivity-weighted radii are the same (see Goad 
et al. 1993), and therefore the measured lag is in- 
dependent of the continuum state. Therefore, if 
after correcting for reverberation effects, the mea- 
sured lag varies with luminosity, then so must its 
effective responsivity. Although Gilbert & Peter- 
son (2003) did not note any evidence for changes 
in rjeff with continuum level, B. Peterson (2003, 
private communication) finds that a detailed anal- 
ysis of the 13 yr campaign on a year-to-year basis 
indicates changes in the line continuum flux slope. 
It is also significant that while BLRs whose lines 
breathe must show an intrinsic Baldwin effect for 
those lines, the converse is not true. For example, 
for the case of 77 = 0.5 everywhere, no breathing 
will occur, but an intrinsic Baldwin effect will be 
present. 

The two triangles in each panel of Figure 3 
mark the positions of the responsivity-weighted 
radii (Rj^) of H/3 and the other three emission lines 
for the two continuum states. In the case of H/3, 
it is larger in the high continuum state by a factor 
of ~ 1.6 relative to the low state. This is about 
the same ratio of luminosity dependent lags ob- 
served in this line in NGC 5548 for the same two 
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continuum states 1.7; Peterson et al. 2002; see 
also Wanders & Peterson 1996; Cackett & Horne 
2005), although we note that the measured lag is 
expected to be smaller than i?,,/c by a factor of a 
few (e.g., see Perez et al. 1992). 

The effective line responsivities for the high and 
low continuum state are given in columns (3) and 
(4) of Table 1. As discussed in §3.1, the low- 
state responsivity is higher in every case, and Fig- 
ure 3 shows that the line's radial responsivity and 
Rri shift to larger distances in higher continuum 
states. For this particular model, the responsivity- 
weighted radius scales as i?,, oc LlonLi with 7 = 
0.15, 0.23, 0.38, and 0.44 for Ha, H/3, He I, and 
He II, respectively. These relations should be 
diagnostic of the physical conditions within the 
BLR. The general relation proposed by Peterson 
et al. (2002), above (7 = 0.5), is expected to 
hold only if the line responsivity is constant over 
a large portion of the BLR with a steep break 
to lower responsivity (see for example. He II in 
Figure 3), and so is not the general case. The 
breathing amplitude (e.g., the value of 7) is ex- 
pected to depend on the responsivity and the ex- 
tent to which the responsivity changes with both 
radius and continuum state in the vicinity of i?,, 
(compare dotted and dashed curves for each of the 
lines in Figure 3). Using the relation between the 
A5100 and A1350 continuum fluxes for NGC 5548 
(Gilbert & Peterson 2003; but see also Peterson 
et al. for a slightly different value), and assuming 
that -L1350 is a proxy for the driving continuum, K. 
Horne (2002, private communication) and Cackett 
& Horne (2005) find a similarly weaker relation be- 
tween the H/3 lag and continuum luminosity state 
for the 13 yr campaign of NGC 5548 as do we. In 
addition, since 7 is related to the responsivity, it 
too changes with continuum level (see also O'Brien 
et al. 1995). 

Finally, wc note that the expected larger line re- 
sponse during low continuum states (§ 3.1) will be 
enhanced by the fact that an effectively smaller 
line-emitting region, due to this breathing, re- 
spond more coherently to a given continuum vari- 
ation. The opposite hold for variations about high 
continuum states. The combined effect can be in- 
vestigated in data sets such as the 13 yr H/3 light 
curve presented in Peterson et al. (2002; their Fig- 
ure 1): more coherent and larger amplitude re- 
sponses in H/3 should generally be found at low 



continuum states, while just the opposite should 
result during high continuum states. One might 
also find differences in line responsivity between 
two comparably low continuum states that imme- 
diately straddle a high continuum state in time, 
due to light travel time effects that are not re- 
moved from the second low state by simply cor- 
recting for the lag. 

3.4. Integrated Line Flux Ratios vs. Con- 
tinuum State 

In this section, we investigate the expected 
changes in the integrated line flux ratios due to 
changes in the continuum state. It should now be 
apparent that such line ratio variations will result 
simply because the responsivities differ between 
the emission lines. Lines with greater responsiv- 
ities vary with greater amplitudes with changing 
continuum state, and so line flux ratios can be ex- 
pected to change as well. Other factors may play 
a role, but responsivity effects will be present re- 
gardless. 

Figure 5 illustrates the Hq/H,/3 line flux ra- 
tio varying across the density-flux plane, from 
roughly 17 (near the triangle at the bottom) to 
approximately 1 (near coordinates: [13.5, 22.5]). 
The classical nebular theory of hydrogen Balmer 
lines (Menzel case B ) does not apply to the gas 
occupying the BLR. For example, at coordinate 
[12.0, 19.5] in the density-flux plane (correspond- 
ing to a distance of 22 It-days in the mean con- 
tinuum state), case B predicts Ha and H/3 fluxes 
that are 20.8 x and 30.4 x larger than that emit- 
ted by the cloud. The case B Ha/H/3 flux ratio 
is 2.35, whereas the ratio predicted by Cloudy is 
3.44. Over a significant portion of the density-fiux 
plane this ratio is inversely proportional to the in- 
cident continuum flux, as expected from the opti- 
cal depth effects mentioned earlier, and in agree- 
ment with the findings of Rees, Netzer, & Fer- 
land (1989; sec also Netzer 1975). This being 
the case, it is easy to imagine this ratio chang- 
ing in response to variations in the incident con- 
tinuum luminosity. Thus, differences in respon- 
sivity among the Balmer lines, as illustrated in 
Figs. 1, 3, 4 and as listed in Table 1, probably 
account for most of the observed steepening in 
the broad-line Balmer decrement as the continuum 
drops into low states, and flattening in high states. 
There have been many such broad line Balmer 
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decrement/continuum state correlations reported 
in the spectra of variable Seyfert 1 galaxies (see, 
for example, Antonucci & Cohen 1983; Ferland, 
Korista, & Peterson 1990; Wamstckcr ct al. 1990; 
Tran, Osterbrock, & Martel 1992). Integrated 
over this model's full broad line-emitting region, 
the Ha/H/? flux ratio changes from 4.9 in the low 
state to 3.7 in the high state. While the partic- 
ulars are model dependent, the trend generally 
is not. The Balmcr decrement's variations with 
the incident continuum level should be sensitive 
to the distribution of clouds in the density-flux 
plane. Furthermore, we suggest that a Seyfert 1 
galaxy's sub-classification as Seyfert 1.0, 1.5, 1.8, 
or 1.9 (as defined in Osterbrock 1989, and ref- 
erences therein) may at least in part be due to 
the relative flux of ionizing photons present within 
their broad-line regions at the time of observa- 
tion. That is, any given spectrum of a Seyfert 
Type 1 AGN may pass through some of these sub- 
classifications as its continuum luminosity fluctu- 
ates from high (Seyfert 1 or 1.5) to low contin- 
uum states (Seyfert 1.8 or 1.9). In high contin- 
uum luminosity states, a larger broad-line contri- 
bution to the total (broad + narrow) Balmer line 
profile will be present and the broad Balmer line 
decrement will be flatter. In low continuum lu- 
minosity states, the broad-line contribution to the 
total Balmer line profile will be relatively much 
smaller, and the broad Balmer line decrement will 
be steeper. NGC 5548 is generally classified as 
Seyfert 1.5. However, the order-of-magnitude con- 
tinuum fluctuations that occurred over the 13 yr 
optical monitoring campaign (Peterson ct al. 2002) 
probably drove the Balmer emission-line spectrum 
from Seyfert type 1.8 in 1992 to perhaps Seyfert 
type 1.0 in 1998. 

As described above, the Balmer lines should 
show less flux variation with continuum state than 
He I, which in turn should show loss variation than 
He II. This trend is generally observed (Peterson & 
Ferland 1986; Ferland, Korista, & Peterson 1990; 
Wamsteker et al. 1990; Dietrich ct al. 1993; Pe- 
terson et al. 1993; KoUatschny 2003). Particu- 
larly interesting is that observed changes in the 
H/3/Hc II flux ratio may, at least in part, due 
to differences in the way these two lines respond 
to the level of the ionizing continuum flux, rather 
than to changes in the shape of the ionizing SED. 
For the model presented here, this ratio drops by 



nearly a factor of 2 from the low to high contin- 
uum states. To a lesser extent, this is likewise true 
for the He II/He I line ratio, which changes by a 
factor of 1.3. Such a responsivity effect may have 
played a role in the observation of a large increase 
in the He II A4686 equivalent width in comparison 
to H/3 during a large increase in the continuum lu- 
minosity in NGC 5548 (Peterson & Ferland 1986) 
and likewise in a similar event reported by Peter- 
son ct al. (2000) in NGC 4051. While the broad 
line He II/H/3 flux ratio is sensitive to the ion- 
izing SED, the usual photon-counting (Zanstra) 
argument used to measure the 13.6 to 54.4 eV 
continuum SED applies strictly when these lines 
are emitted in gas whose properties allow for re- 
sponsivitics rj^ff « 1, as would be true if case B 
applied for both lines. However, these conditions 
do not apply for the Balmer lines emitted within 
the BLR, but they can be found within the narrow 
line regions of AGNs (see the EW distributions of 
the recombination lines in Ferguson et al. 1997) or 
within galactic nebulae. On the other hand, the 
He II emissivity does not stray too far from that 
expected from case B for most conditions expected 
in the BLR, as indicated by its comparatively fea- 
tureless EW distribution in the density-flux plane 
for log$H < 21 (see also Bottorff et al. 2002). 

Finally, as a matter of general interest we note 
that the Lya/H/3 flux ratio should correlate with 
the continuum flux, and apparently this is the case 
in the well-studied AGN NGC 5548 (Wamsteker 
et al. 1990). The EW contours of Lya in the 
density- flux plane (sec Figure 1, bottom right) do 
not show as strong continuum flux dependencies as 
do the Balmer lines, and so ry(Lya) > 7?(H/3) can 
be expected, with larger differences to come from 
gas nearer to the continuum source where the re- 
sponsivity of H/3 is small (see also Figure 4) . The 
combined conclusions of Pogge & Peterson (1992) 
and Gilbert & Peterson (2003) would indicate an 
observed relation of L(Lya)/L(H/3) oc L?35o"°'^^ 
in NGC 5548. Our simple model (see Table 1) 
predicts a similar relation _L(Lya:) /L(Rf3) cx £1350, 
which is steeper in high continuum states and shal- 
lower in low ones. Shields & Ferland (1993) dis- 
cussed some of the mechanisms behind the de- 
struction of Lya emitted by broad-line clouds. 
Pogge & Peterson also found that in Fairall 9 
rjeffiLya) « 0.71 for a factor of ^10 change in the 
observed A1338 continuum, similar to the value re- 
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ported here (0.74, Table 1) for our simple model. 
Figure 4 shows that r7(Lya) is nearly independent 
of gas density and incident photon flux over a wide 
range in these parameters. 

3.5. Line Profiles and Their Variations 
with Continuum State 

While the dynamics of the gas emitting the 
broad emission lines are not yet understood, the 
importance of a virial component seems to be 
reasonably well established (Peterson & Wandel 
2000). Here, we have adopted the simple relation- 
ship 

v{r) = {GM/ry/^, (3) 

assuming randomly inclined circular orbits and a 
central mass of 8.3 x 10'' solar masses, chosen to 
approximately reproduce the observed broad-line 
FWHM of H/3 in NGC 5548. This choice of dy- 
namical model is motivated by simplicity for pur- 
poses of demonstration. 

The low-state contimium emission-line profiles 
of the four lines under study, as predicted by 
our simple LOG model line emissivity and ve- 
locity field, are shown in Figure 6, normalized 
to their peak intensities. We note that the flat 
tops to these profiles are artifacts of the assump- 
tion of circular orbits, and their extents are de- 
termined by the local Keplerian velocity at the 
distance of the model's outer radius (140 It-days) 
given the assumed central mass. The EW distri- 
butions shown in Figure 1, and most notably the 
continuum flux at the onset of significant ther- 
malization, are roughly reflected in the widths of 
the emission-line proflles in Figure 6, as well as 
in their observed lags (Dietrich et al. 1993; Pe- 
terson & Wandel 1999, 2000; Fromerth & Meha 
2000). He II is the broadest line (with the short- 
est lag), followed by He I, then H/3, and finally 
Ha (with the longest lag). This is as generally 
observed in broad-line AGNs (Stirpe 1990, 1991; 
Grenshaw 1986, Shuder 1982, 1984; Osterbrock & 
Shuder 1982; Boroson & Green 1992; KoUatschny 
2003). Krolik et al. (1991), Peterson & Wan- 
del (1999, 2000), and Onken & Peterson (2002) 
present highly correlated lag line width relations 
that indicate virial-like motions in the vicinity of a 
supermassive central object. The FWHM of these 
four emission lines, as well as H7 and Lya, are 
listed in columns (5) and (6) of Table 1 for the 



high and low continuum states, respectively. The 
trend in line widths is in general agreement with 
those found in the models presented by Rees, Net- 
zer, & Ferland (1989). 

Since the lines under study here are recombi- 
nation lines involving just one or two electrons, 
they are emissive over very broad areas within the 
density-flux plane, unlike the collisionally excited 
metal lines which are very sensitive to the ion- 
ization parameter. The Balmer and He I lines 
span essentially equal areas of signiflcant emissiv- 
ity within the density-flux plane without regard to 
ionization parameter, and that of He II does like- 
wise for logUn ^ —3.5 (see Figure 1). While it is 
also true that a cloud of a given column density 
will become fully Hc++ at a higher ionization pa- 
rameter than for H+, this difference is small (see 
Figure 1). By far most of the differences in emis- 
sivity among the Balmer and He I lines arc due to 
the optical depth effects discussed above, and this 
is also true for the He II lines except for clouds 
with logUn ^ —3.5. Stated another way, if all of 
these lines were emitted with their case B emissiv- 
ities within the BLR, much smaller differences in 
their profiles and lags would be expected. 

Figure 7 shows the long term (global) respon- 
sivities for each of the four emission lines as a func- 
tion of radial velocity across their line profiles, for 
our simple velocity field. Because of the larger 
response in the line cores, the broad emission- 
line component of the Balmer lines should become 
narrower/peakier in higher continuum states, and 
broader/less peaky during low continuum states 
(cf.. Ha variations in NGG 5548; Stirpe et al. 
1988); refer also to the model high and low con- 
tinuum state FWHMs as listed in columns (5) and 
(6) in Table 1. In fact for this particular model, 
a continuum state that is a factor of 12 x brighter 
than the low state'' produces a FWHM in H/3 that 
is 1000 kms"-"^ narrower than that in the low state. 
The response across the profile of He II is much 
higher, although again this emission line would be 
expected to be somewhat narrower during higher 
continuum states. The response across the pro- 
file of He I is not as broad as He II and not as 
peaky as the Balmer lines. In every case, we ex- 

'Such a high state in NGC 5548 was recorded in 1984 
(Wamsteker et al. 1990) and again in 1998 (Peterson et al. 
2002). 
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pect the response in the line cores to be stronger 
than in the line wings. The relative change in 
the FWHM of the emission Hne with continuum 
state is effectively a reflection of the "breathing" 
BLR as caused by changes in the line respon- 
sivity, as we discussed in § 3.3. O'Brien, Goad, 
& Gondhalckar (1995) similarly linked non-linear 
emission-line response to line profile variability, al- 
though this general conclusion was based on single 
power law cloud models. 

3.6. Radial Velocity-Dependent Emission- 
Line Ratios and Their Variations with 
Continuum State 

Since the responsivity of these lines is expected 
to be radially dependent and this function dif- 
fers between emission-line species, we expect ra- 
dial velocity-dependent variations in the line in- 
tensity ratios. These effects are illustrated in Fig- 
ure 8, which shows the velocity-dependent Balmer 
decrement (Ha/H/3) for the low-state {dotted line) 
and high-state {solid line) continuum. Note that 
in both cases the Balmer decrement is steeper in 
the line core than in the line wings, and the Balmer 
decrement shows larger changes with continuum 
level in the line core (from 5.6 to 4.0), than in the 
line wings (from 3.0 to 2.7). The first of these ef- 
fects has been well-established (see, for example, 
Shuder 1982, 1984; Osterbrock & Shuder 1982; 
van Groningen 1984, 1987; Crenshaw 1986; Stirpe 
1990, 1991; Kollatschny 2003), and we are aware 
of just one study of the time-dependent variabil- 
ity of the Balmer decrement as a function of ra- 
dial velocity (Stirpe, de Bruyn, & van Groningen 
1988). Both behaviors imply that the velocity field 
must decrease with increasing radius from the cen- 
tral source. As just presented above, a weaker re- 
sponse in an individual line's wings in comparison 
to its core (Fig. 7) would also argue for such a 
velocity field. These findings are consistent with 
recent reverberation studies that find a relation- 
ship between the emission-line lag and rms profile 
width for multiple lines in individual AGNs, which 
implies that the velocity field is dominated by the 
central massive object (Peterson & Wandel 1999, 
2000; Fromerth & Melia 2000). 

The results presented here also predict a smaller 
H/3/He I intensity ratio in the line wings than in 
the core (Fig. 8). The EW contours in Figure 1 
coupled with a velocity field that diminishes with 



increasing distance from the central source provide 
possible explanations for this generally observed 
phenomenon (Shuder 1982, 1984; Osterbrock & 
Shuder 1982; Crenshaw 1986; van Groningen 1984, 
1987; Rees, Netzer, & Ferland 1989; Kollatschny 
2003). These lines' EWs peak up at similarly high 
densities, but the He I line peaks up at signif- 
icantly higher photon fiuxes (smaller radii), be- 
cause of its lower optical depth. This is responsi- 
ble for this line's relatively stronger wings (see also 
Netzer 1978; Kwan 1984; Carroll 1985; Carroll & 
Kwan 1985). 

Whether or not the specific dependencies of the 
continuum state/responsivity effects described in 
this section are present in the spectra of NGC 5548 
or other variable Seyfert 1 AGNs is not par- 
ticularly important. We discuss the sensitivity 
of line responsivity to model assumptions in the 
Appendix. Instead, the emphasis of this paper 
lies mainly in illuminating the mostly untapped 
wealth of the potentially powerful physical diag- 
nostics present in high quality spectrophotometric 
data sets of the optical recombination lines, and 
that their emissivities' sensitivity to the contin- 
uum flux makes them particularly effective diag- 
nostics of the radially dependent physical condi- 
tions within the BLR. 

4. DISCUSSION 

4.1. Optically Thin Gas in the So-Called 
Very Broad Line Region 

That the BLR may contain a significant com- 
ponent of ionized gas that is optically thin to 
hydrogen-ionizing photons has been suggested by 
investigators over many years. In nearly every 
case there is an association of this gas with that 
emitting the very broad wings of the lines. Most 
recently, Sulentic et al. (2000) reported that in 
the quasar PG 1416-129 the FWHM of the broad 
emission-line profile of H/3 increased by 50% from a 
high to a low continuum state separated in time by 
10 yr. Their interpretation of this is that the wings 
of the Balmer lines are dominated by emission 
from an optically thin "very broad line region" 
(VBLR), while the core of the line is dominated 
by more highly variable optically thick gas. These 
are usually thought of as two distinct line-emitting 
regions. Ferland, Korista, & Peterson (1990), Pe- 
terson et al. (1993), and Corbin & Smith (2000) 
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made similar arguments based on B aimer line - 
continuum variations, while Gondhalekar (1987, 
1990), O'Brien, Zheng, & Wilson (1989), & Perez, 
Pcnston, & Moles (1989) foimd less variable wings 
in Lya and C IV A 1549 of higher redshift quasars. 
Kassebaum et al. (1997) found similarly stronger 
variations in the line core than wings of H/3 in 
Mrk 335. It should be noted that while the con- 
clusions regarding the core versus wing variations 
from monitoring campaign studies such as Peter- 
son et al. (1993) and Kassebaum et al. (1997) were 
based on analyses of time-resolved spectral obser- 
vations of the continuum and H/3, most of the 
above studies relied on very few, often just two, 
observation epochs. Without sufficient time reso- 
lution of the variability, it is difficult to assess the 
actual line response. Additional arguments for an 
optically thin line-emitting region have been put 
forth through the analyses of line flux ratios in 
the profile wings: optical Balmer and He II lines 
(Marziani & Sulentic 1993; Corbin 1997), Balmer 
and C IV A1549 (Corbin 1995), Balmer and Lya 
line profiles (Zheng 1992), and Balmer and O I 
A8446 profiles (Morris & Ward 1989). 

While at any given continuum state the broad 
emission-line region may well contain clouds fully 
ionized in hydrogen, as the present model does, 
the results presented here suggest that there is 
no great need for a major line-emitting compo- 
nent consisting of vast amounts of optically thin 
gas. Most of the effects mentioned just above 
regarding the Balmer and helium recombination 
lines can be explained by the eff'ects we described 
in § 3: weak Balmer line emission and response 
need not originate in clouds optically thin to the 
Lyman continuum. While an emission line such 
as Ne VIII A774 is likely to be emitted mainly 
within clouds that lack a hydrogen ionization front 
(e.g., see Hamann et al. 1998; Korista et al. 1997), 
we have also pointed out that by their very na- 
ture fully-ionized clouds are inefficient emitters of 
hydrogen recombination lines (refer again to Fig- 
ure 1), and so arc unlikely to important en- 
ergetically even if their covering fraction is large. 
This is especially true in the face of a variable 
ionizing continuum source; clouds will more likely 
be either ionization-bounded or too overionized to 
emit Balmer lines of any significance. It is also 
significant that the matter-bounded clouds in Fig- 



ure 4 either have small positive (lightest grayY or 
negative (white) responsivities. Even if present, 
the latter type clouds never make their presence 
known in these emission lines, which is not surpris- 
ing given their extremely low emitting efficiencies. 

In regards to the comparisons of Balmer and 
Lya line profiles, as mentioned in § 3.4, the EW 
contours of Lya in the density-flux plane (see Fig- 
ure 1, bottom right) do not show the strong con- 
tinuum flux dependencies of the Balmer fines. We 
suspect that this, along with the presence of a ra- 
dially decreasing velocity field, accounts for the 
very large Lya/H/3 ratios found in the emission- 
line wings by Zheng (1992), rather than these ra- 
tios being due to a major constituent of optically 
thin emitting gas. 

Finally, we also point out that many of the 
above studies assumed that case B emissivities ap- 
ply to the hydrogen lines emitted by optically thick 
clouds in the BLR and attributed significant de- 
viations from expected ratios to an optically thin 
emitting region. We conclude that a separate op- 
tically thin VBLR may be unnecessary, and in any 
case for the Balmer lines is probably energetically 
unimportant. The wings of the hydrogen recom- 
bination lines likely represent nothing more than 
the inner BLR where the Balmer line emission is 
inefficient and the responsivity low. 

4.2. Effects of an SED that Varies with 
Continuum Level 

Paltani & Courvoisier (1994) and Paltani & 
Walter (1996) showed unambiguously that the UV 
continuum becomes bluer when the flux is higher 
in variable AGNs. This effect has also been found 
in the Sloan Digital Sky Survey quasar sample 
(Vanden Berk et al. 2004). In the particular 
case of NGC 5548, Maoz ct al. (1993) reported 
that, after accounting for the Balmer continuum 
and UV Fe II emission, the 1800-2400 A con- 
tinuum becomes bluer with increased continuum 
flux. Additionally, Peterson et al. (2002) and 
Gilbert & Peterson (2003) found that after care- 
ful removal of the contaminating non- variable stel- 
lar optical continuum flux, the UV-optical contin- 
uum in NGC 5548 becomes bluer as it becomes 



The only clouds in Figure 4 whose hydrogen line responsiv- 
ities are lie above the diagonal line log(f/Hc) ^ 10.25 
and also within the lightest gray shaded region. 
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brighter: F5100 oc -Fj'sfo- Extrapolation of this 
relation to higher energies might indicate a hard- 
ening of the incident ionizing continuum during 
brighter contimmm states. Marshall ct al. (1997) 
found larger amplitude variations in the extreme 
ultraviolet (~ 100 A) than at 1350 A during a 
week-long period of the 1993 NGC 5548 monitor- 
ing campaign. Alternatively, as suggested by Pal- 
tani & Walter (1996), it may indicate the presence 
of a constant (or more slowly varying) continuum 
or pseudocontinuum component whose contribu- 
tion is greater at longer wavelengths. Korista & 
Goad (2001) found that 20%-30% of the observed 
effect may be due to reverberation of the diffuse 
continuum emission of the broad line clouds. In 
the simulations presented here, the SED of the 
continuum incident on the broad-line clouds is the 
same for both high and low continuum states. How 
an emission line's responsivity would change with 
an incident SED that changed shape as the UV- 
optical continuum level varied will depend on the 
details of how the flux of photons important to 
the creation and destruction of the line varies. We 
briefly discuss these effects here. 

The shape of the ionizing continuum deter- 
mines the heating for given values of gas den- 
sity, ionizing photon flux, and chemical abun- 
dances. All else being equal, harder incident con- 
tinua result in higher electron temperatures, and 
so major cooling lines (e.g., C IV A1549) might 
be expected to increase their energy output in 
response. One should also keep in mind that 
the incident flux at energies corresponding to the 
Balmer continuum is important to the destruction 
of hydrogen emission lines for the high gas densi- 
ties encountered in the BLR (Shields & Ferland 
1993). Also, harder continua during brighter UV- 
opt continuum states might also mean compara- 
tively more photons above 1 ryd and/or 4 ryd, 
favoring the production of hydrogen and helium 
recombination lines. Kaspi & Netzer (1999) con- 
sidered an SED whose break near 4 ryd moves to 
higher energies for brighter UV continuum states 
(their Figure 8d), with the effect of increasing 
the 4 ryd/1 ryd photon flux ratio with increas- 
ing A1350 continuum luminosity. However, in 
this variable SED scheme, the flux above 20 ryd 
(270 eV) remained constant. A comparison of 
their Figs. 7 (fixed SED) and 9 (variable SED) that 
show their model predictions of five UV emission 



lines' light curves versus observations illustrates 
some of the aforementioned effects of a variable 
continuum SED on line responsivity. While they 
included He II A 1640 in their simulations, such was 
not the case for the Balmer lines. Their prescrip- 
tion for a variable continuum SED increased the 
responsivities of the He II and Lya lines, whereas 
those of C IV A1549, C III] A1909, and Mg II A2800 
decreased, presumably because the variations in 
the continuum beyond ^150 cV were smaller in 
their variable SED scheme than for the constant 
SED. It is also apparent that the lags of He II 
and C IV increased in the higher continuum states, 
because of a shift to larger responsivity- weighted 
radii. 

4.3. The Importance of the Recombina- 
tion Lines to the Method of Quasar 
Tomography 

The method of "quasar tomography," proposed 
by Horne et al. (2003), is a unification of re- 
verberation mapping and photoionization physics. 
Its goal is to place physical constraints on the 
recovered two-dimensional transfer function map 
^'(t, t;) of the broad emission lines, as well as to 
recover the gas distribution described by a five- 
dimensional map f{r,9, uh, Nh, v), where r is the 
time delay, v is the observed radial velocity, r 
is the distance from the continuum source, and 
9 is the angle from the line of sight. The sen- 
sitivity of the emissivities (EW contours in Fig- 
ure 1) of the optical recombination lines to the 
continuum flux, and the resulting consequences to 
their responsivities, make these lines important 
candidates for inclusion in the quasar tomogra- 
phy. Horne et al. found that when applied to the 
mostly coUisionally excited UV emission lines, a 
certain amount of ambiguity arises in the recovery 
of the maps due to the weak dependencies in their 
EWs (and so responsivities) along lines of constant 
ionization parameter (cx ^H/nn) corresponding 
to their maximum emissivities. Their emissivities 
and responsivities are largely dependent on Uh, 
with smaller dependencies in gas density due to 
de-excitation or thermalization. The optical re- 
combination lines, along with Mg II A2800, offer 
additional constraints in determining the physical 
conditions within the BLR, in that their emissivi- 
ties and responsivities are dependent on the local 
continuum flux and so radius for a fixed central 
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continuum luminosity, as well as changes in this 
flux as the central continuum source varies. 

5. SUMMARY 

The response of the optical recombination lines 
to continuum variations, especially the Balmer 
lines, depends critically on the time-averaged pho- 
ton flux incident on a spatially extended BLR. 
This strong dependence exists because of the en- 
hanced excited state populations (and so optical 
depths) that occur at higher continuum flux levels, 
given the high gas densities found in the BLRs of 
AGNs. The emissivities, and hence responsivities, 
1], of these linos arc generally strongly anticorre- 
lated with the continuum flux, so r] changes with 
distance within the BLR for a fixed continuum lu- 
minosity and with time as the continuum source 
varies. Below are the major implications, most of 
which provide explanations for observed relation- 
ships: 

1. We predict an intrinsic Baldwin effect due to 

this non- linear response, specifically r) < I, 
in most emission lines, including the Balmer 
lines, even after accounting for geometric di- 
lution duo to reverberation. This is a finding 
that has only recently been observationally 
confirmed. Regarding the optical recombi- 
nation lines, this effect should be small for 
He II A4686 (with responsivity r; « 1), but 
significant for the Balmer lines. 

2. We predict that the higlica- order Balmer 
lines arc more responsive. rj{H"/) > r](H(3) > 
r]{Ha), as is generally observed, and so the 
broad-line Balmer decrement will steepen in 
low continuum states, and flatten in high 
states. Moreover, He II A4686 should be 
substantially more responsive (0.9 ^ f] < 1) 
than the Balmer lines, and so large varia- 
tions in the broad He II/H/3 integrated flux 
ratio that are correlated with the continuum 
state can be expected, even in the absence 
of a varying continuum SED. 

3. A Seyfert 1 galaxy's sub-classiflcation as 
Seyfert 1.0, 1.5, 1.8, or 1.9 may at least 
in part be due to the relative flux of ioniz- 
ing photons present within their broad-line 
regions at the time of observation. Like- 
wise, any given spectrum of a Seyfert Type 1 



AGN may pass through some of these s\ib- 
classifications as its continuum luminosity 
fluctuates from high (Seyfert 1 or 1.5) to low 
continuum states (Seyfert 1.8 or 1.9). 

4. The emission-line responsivity, ry, places an 
additional important constraint on the run 
of physical conditions within the BLR (e.g., 
gas density, column density, covering frac- 
tion) as a function of radius over and above 
that of the line emissivity and lag. This is 
because for the optical recombination lines, 
the mean responsivity declines with increas- 
ing incident photon flux. Thus, for a partic- 
ular continuum luminosity state, a smaller 
emissivity-weighted radius is accompanied 
by a lower effective responsivity. 

5. Greater responsivity in the optical recom- 
bination lines is expected during low con- 
tinuum states, and lower responsivity is ex- 
pected during high continuum states. The 
characteristic size of the BLR as measured 
by the emission-line lag also will track the 
continuum luminosity level, a phenomenon 
known as "breathing." Responsivities that 
are functions of the continuum level and 
hence time result in non-stationary line 
transfer functions. In this case the time 
dependence is due to changes in the time- 
averaged photon flux incident on a spatially 
extended BLR. The stronger response due 
to an elevated responsivity at low contin- 
uum levels will be enhanced by the greater 
coherence of an effectively smaller BLR due 
to breathing. This time-dependence in the 
line transfer functions should also manifest 
itself in changes in the lino profile widths, 
with the broad-line component becoming 
narrower /peakier in high continuum states. 

6. Because the emissivities and responsivities 
of the optical recombination linos are anti- 
correlated with the incident continuum flux, 
observations that Hq:/H/3 and H/3/He I flux 
ratios and the Balmer line responsivity de- 
crease from the core to the wings indicate 
that the BLR velocity field diminishes with 
increasing distance from the central contin- 
uum source. These findings are consistent 
with recent reverberation studies that find a 
strong correlation between the emission-line 
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lag and rms profile width for multiple lines 
in individual AGNs, which implies that the 
gas kinematics are dominated by the central 
massive object. 

7. Larger variations in the Balmer decrement 
with continuum state will be found in the 
broad-line cores rather than in the wings. 
Observations of such would also indicate 
that the BLR velocity field diminishes with 
increasing distance from the central contin- 
uum source. 

8. Differences in broad-line widths (whether 
rms or mean profile) and lags among the 
Balmer and helium recombination emission 
lines are in large part due to a radial (contin- 
uum fiux-dependent) stratification of optical 
depth effects in these lines. 

9. Much of the observed behavior of the opti- 
cal recombination lines (e.g., the weak re- 
sponse of the Balmer line wings) previously 
attributed to a substantial contribution from 
matter-bounded clouds at small BLR radii 
(i.e., the so-called VBLR) may be explained 
by the photon flux dependencies of these 
lines' emissivities and responsivities within 
ionization-bounded clouds. 

These findings should be general, largely inde- 
pendent of any specifically adopted model of the 
BLR. However, the details of the optical recom- 
bination lines' responsivities will be sensitive to 
and so diagnostic of specific physical conditions 
within the BLRs of AGNs. In addition to the 
recombination lines investigated here, we suggest 
that studies of the IR Paschen and Bracket lines 
should prove fruitful in further constraining the 
BLR, particularly in its outer regions. 
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Appendix: SENSITIVITY OF EMISSION LINE RESPONSIVITY TO MODEL ASSUMPTIONS 
We briefly discuss the point that most of the major trends in hne responsivity reported here should not 
depend critically on the particular choice of model for the broad emission-line region. As discussed in § 2, 
most of these effects can be understood from Figure 1, which shows quite generally the dependence of an 
emission line's emissivity on cloud density and incident ionizing photon flux, irrespective of how these clouds 
(or emitting entities) are distributed or other details. 

A. Cloud Column Densities 

Korista et al. (1997) discusses and illustrates the effects of cloud column density on the line emissivities, 
and thus EW contours in Figure 1. The present simulations assume a constant cloud column density of 
10^^ cm~^. In general, the computed integrated AGN line spectrum is relatively insensitive to this parameter 
for column densities exceeding 10^^ cm~^ (see also Goad & Koratkar 1998), except for the high ionization 
lines of O VI A1034 and Ne VIII A774, which require colmnn densities > 10^^-^ cm"^ and > 10^^ cm"^ to fully 
form within the cloud (Korista et al. 1997; Hamann et al. 1998). To the extent that cloud column density 
has a significant effect on the luminosities or responsivities of the emission lines plotted in Figure 1, the 
clouds must lie in the vicinity of the central diagonal boundary representing the steep downward gradient in 
EW (i.e., near log $i/ — log riij ~ 10.25). These would be the clouds affected by the passage of an ionization 
front. If the preponderance of clouds have a smaller column density than those computed here, this diagonal 
boundary would shift toward the lower right corners of the density- flux planes in Figure 1. For higher column 
density clouds, this boundary would shift toward the upper left. The emissivities of clouds lying far from this 
boundary in the density- flux plane would be essentially unaffected. This is clearly illustrated in Figure 3 of 
Korista et al. (1997). From Figure 4 and Eq. 2 in the present paper, it can therefore be seen that changing 
the cloud column density by an order of magnitude in either direction from the adopted value will not 
change the overall trends of the EW and responsivity of the ionization-bounded clouds being proportional 
to the incident continuum flux for the four emission lines studied here. Because of their inefficient natures, 
matter-bounded clouds cannot contribute substantially to the hydrogen recombination line luminosity (or 
responsivity), unless they have much larger covering fractions than the optically thick clouds at the same 
radius. Even then, their importance is likely to be minor, except perhaps in the extreme wings given a 
velocity field that declines with distance from the central source. 

B. Incident Continuum SED, Gas Abundances 

Different assumptions about the incident continuum SED or gas abundances do not change the line EW 
distributions in Figure 1 very much (see Korista et al. 1997), and so would have little effect on these lines' 
responsivities. 

C. Optical Depth and Radiative Transfer Effects 

If the optical depths in the Balmer lines were significantly reduced by the presence of large extra-thermal 
velocities that might occur in an accretion disk wind (e.g., Murray & Chiang 1997), the result would be the 
EW contours in Figure 1 peaking up at higher incident fluxes and densities and the Balmer line-emitting 
region being more strongly influenced by gas with high responsivity (?7 ~ 1^2). Since the Balmer lines have 
reported responsivities less than 1 even after correcting for geometric dilution, this would constrain models 
with large extra-thermal velocities and/or Balmer emission from large outer radii. Other consequences of 
signiflcant extra-thermal velocities might be smaller predicted lags and broader predicted line profiles (for a 
fixed central mass and incident continuum luminosity). If in general the radiative transfer of the Balmer lines 
is currently less than adequately handled (as may be the case; see Netzer et al. 1995), the optical depth of 
Ha should nevertheless remain greater than that of H/J, and those of the helium lines should remain smaller 
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than the Balmcr hncs. Thus, the relative senses of the EW contours of the Balmer hues in Figure 1 should 
be preserved. The sense of the Balmer decrement variability as a function of continuum state should also be 
preserved, as should the Ra/Hp intensity ratio across the line profile: provided that the bulk of the velocity 
field diminishes with increasing distance from the central source. 

D. Possible Alternate Trends in Lag and Line Response 

Most of the extra rcsponsivity in H/9 over Ha in Figure 7 lies within the line core whose contributions 
are dominated by gas lying (and responding) at larger distances. As a consequence, the difference in the 
measured lags of these two lines may not be as large as their EW distributions in Figure 1 might indicate. 
That is, their responsivity-wcightcd radii may be closer in size than their cmissivity-wcighted radii (Goad 
et al. 1993). In fact, depending on the actual emissivity and responsivity distributions of the two lines, it is 
possible for H/3 to lag slightly behind Ha. In other words the total light profile may be broader in H/3 than 
Ha and the other way around in the rms variation profile, with Ha having a slightly smaller lag. The small 
differences in the Balmer lines' beaming functions (O'Brien et al. 1994) may also play a role in altering the 
relative responsivity-weighted radii, and thus relative lags of these emission lines. 
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Log n(H) Log n(H) 

Fig. 1. — Contours of logT40\ (EW) for six emission lines, referenced to the incident continuum at 1215 A, 
are shown as a function of the hydrogen number density and flux of hydrogen-ionizing photons for full source 
coverage at every position. The total hydrogen column density is 10^'^ cm~^. The EW is in direct proportion 
to the continuum reprocessing efficiency. The smallest decade contour (generally rimning diagonally throiigh 
the centers of the diagrams) corresponds to 1 A, each solid line is f decade, and dash-dotted lines represent 
0.1 decade steps. The contours decrease monotonically from the peak {triangle) to the 1 A contour, and the 
blank regions in the upper left portion of the panels have EWs less than 1 A. The star in each panel is a 
reference point marking the old "standard BLR" parameters (Davidson & Netzer 1979). 
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Fig. 2. — Effective logarithmic H/3 surface flux (ergs s^^cm^^) versus the logarithm of the distance from the 
ionizing continuum source for the two continuum states. The H/J emissivities of clouds with a range of gas 
densities have been summed at each radius. 
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Fig. 3. — Radial responsivities of the four recombination lines under study. The solid line represents the 
responsivity for the factor of 8.2 change in continuum flux. The dotted line represents the same for 10% 
variations about the low continuum state, and the dashed line for 10% variations about the high continuum 
state. The responsivity-weighted radius is marked with a triangle for each of the latter two cases. 
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Fig. 4. — Same as Fig. 1, but with a superposed gray scale in responsivity. White represents rj < 0. The 
lightest shade of gray represents Q < rj < 0.5, medium gray 0.5 < rj < 1.0, medium dark gray 1.0 < < 1.5, 
dark gray 1.5 < 77 < 2.0 (appearing in He II A4686, H/3, H7 only), and black r] > 2 (appearing only in the 
panel for He II A4686, which shows all of the shades). These responsivities were computed for 0.45 decades 
(factor of 2.8) variations in the continuum luminosity about the adopted mean continuum state. 
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Log n(H) 

Fig. 5. — Logarithmic flux ratios of Ha/H/3 in the density-flux plane. The sohd contour near the bottom of 
the diagram represents a flux ratio of 10 (1.0 in the log), the dash-dotted contours are 0.1 decade intervals, 
and the solid contour in the upper right corner represents a flux ratio of 1.0 (0.0 in the log). The triangle 
lies at the peak ratio of w 17. The star is a reference point marking the old "standard BLR" parameters of 
Davidson & Netzer (1979), and has a value of w 4.2. The contours end abruptly at the central diagonal solid 
line, where the equivalent widths of the two lines become very small due to the extremely high ionization 
state of the gas. 
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Fig. 7. — Model emission-line responsivities presented as functions of radial velocity for the four recombina- 
tion lines under study (§ 3.5). 
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Fig. 8. — Ha/H/3 and H/3/He I A5876 fliix ratios presented as functions of radial velocity in the low and high 
continuum states (§ 3.6). In both cases, the smaller ratios are found in the high continuum states because 
of the stronger response of the line in the ratio's denominator. 
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Table 1: Model Effective Line Responsivity & Profile Variability 
Emission Line 7?e// (global) r;e//(hi) 7?e//(lo) FWHM(hi) (km/s) FWHM(lo) (km/s) 



Ha A6563 0.52 0.43 0.63 4200 4600 

H/3A4861 0.64 0.54 0.77 4510 5230 

H7 A4340 0.68 0.58 0.79 4570 5250 

He I A5876 0.83 0.77 0.87 5780 7350 

HeIlA4686 0.92 0.87 0.97 7590 9480 

Lya A1216 0.74 0.71 0.77 5180 5280 
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